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ABSTRACT 


\ 

\ 

l 

The  IBM  Company's  Continuous  System  Modeling  Program 
was  used  to  simulate  the  lateral  and  directional  flight 
control  systems  of  the  F/A-18  aircraft.  The  model  is 
designed  for  use  in  studies  of  high  angle-of -attack  maneuver 
ing  flight  and  is  restricted  to  the  Auto  Flaps  Up  mode  of 
operation.  The  model  accepts  simulated  pilot  stick  and 
rudder  inputs,  air  data  information,  and  rate  gyro,  angle- 
of-attack,  and  acceleration  feedback  signals.  Outputs  are 
differential  stabilator,  differential  leading-edge  and 
trailing  edge  flap,  aileron,  and  rudder  deflections. 

Typical  in.put  values  are  used  to  validate  the  model, 
generating  output  control  surface  deflections  which  corres¬ 
pond  to  those  expected  for  the  F/A-18  aircraft. 
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I.  INTRODUCTION 


In  this  thesis,  work  completed  at  the  Department  of 
Aeronautics  of  the  Naval  Postgraduate  School  as  part  of  that 
Department's  research  into  the  development  of  advanced 
control  concepts  using  digital  electronic  flight  control  in 
U.  S.  Navy  tactical  combat  aircraft  is  described.  In  con¬ 
current  work.  Carter  [Ref.  1],  discusses  in  part  the 
motivation  for  and  the  scope  of  this  research  program. 
Briefly  reviewing  Reference  1,  it  is  desired  to  expand  the 
application  of  digital  flight  control  technology  to  the 
following  areas: 

1.  Investigation  of  active  control  prevention  of  departure 
from  controlled  flight. 

2.  Comparative  testing  of  new  control  law  algorithms. 

3.  Evaluation  of  modern  control  techniques  such  as  optimal 
control,  observers,  and  model  following  control. 

4.  Simulation  of  degraded  flight  conditions  for  combat 
survivability  studies. 

As  discussed  in  Reference  1,  computer  simulation  of  a 
digital  fly-by-wire  aircraft  was  deemed  the  best  initial 
approach  in  studying  the  topics  listed  above.  The  Navy/ 
McDonnell  Douglas  F/A-18  is  such  an  aircraft,  and  in  fact  is 
the  first  U.  S.  aircraft  in  production  to  use  a  digital 
fly-by-wire  control  system.  In  Reference  1,  the  computer 
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modeling  of  the  longitudinal  axis  control  system  of  the 
F/A-18  is  described.  In  this  thesis,  the  model  is  extended 
to  the  lateral-directional  axes  to  complete  the  control 
system  model.  Future  incorporation  of  an  F/A-18  aerodynamics 
simulation  by  Raithel  (Ref.  2]  and  non-linear  equations  of 
motion  will  result  in  a  complete  "flying"  computer  model  of 
the  F/A-18. 

The  computer  model  was  developed  on  the  Naval  Postgrad¬ 
uate  School's  IBM  Model  3033  general  purpose  mainframe 
computer.  The  code  is  written  for  the  IBM  Company’s 
Continuous  System  Modeling  Program  (CSMP)  (Ref.  3],  which  is 
a  Fortran  application  program  designed  to  simulate  dynamic 
systems.  The  text,  "A  Guide  To  Using  CSMP"- by  Speckhart  and 
Green  [Ref.  4],  provides  most  of  the  documentation  necessary 
to  use  CSMP.  CSMP  has  several  advantages  over  a  conventional 
Fortran  program  for  this  simulation.  First,  CSMP  includes 
34  built-in  functions  which  serve  to  model  most  of  the  block 
diagrams  encountered  in  control  systems  engineering.  These 
functions  are  analogous  to  Fortran  functions.  Second,  CSMP 
can  accommodate  user-defined  functions,  called  macros,  and 
can  also  call  standard  Fortran  subroutines  and  functions. 
Third,  output  formatting,  either  printed  or  graphical,  is 
handled  by  CSMP  using  just  a  few  statements.  Finally,  the 
time  base  and  numerical  integration  routines  are  provided 
by  CSMP.  One  might  consider  the  primary  disadvantage  of 
CSMP  to  be  lack  of  user  control  over  the  built-in  functions 


and  numerical  methods.  For  this  simulation  CSMP  proved 
to  be  a  very  effective  tool. 

In  Chapter  II,  an  overview  of  the  F/A-18  flight  control 
system  is  provided  and  the  assumptions  made  and  limitations 
of  the  modeling  process  are  discussed.  In  Chapter  III,  the 
methodology  and  nomenclature  of  the  computer  program  are 
outlined.  In  Chapter  IV,  the  tests  and  results  used  to 
validate  the  model  are  presented.  Conclusions  and  recommen¬ 
dations  are  in  Chapter  V.  Appendix  A  contains  block  diagrams 
of  the  F/A-18  flight  control  system  as  modeled,  including  the 
previous  longitudinal  model  developed  by  Carter  [1],  Changes 
in  nomenclature  and  arrangement  to  Carter's  model  have  been 
made  in  incorporating  it  into  the  current  model;  these  are 
essentially  cosmetic.  Appendix  B  is  the  computer  program 
listing  in  CSMP  of  the  F/A-18  flight  control  system  as 
modeled,  again,  readers  of  Carter's  program  will  note  some 
changes  in  its  form  here  to  suit  the  overall  program. 
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II .  F/A-18  FLIGHT  CONTROL  SYSTEM 

The  F/A-18  flight  control  system  is  described  in  detail 
in  the  McDonnell  Aircraft  Company's  F/A-18A  Flight  Control 
System  Design  Report  [Ref.  5].  In  this  chapter,  basics  of 
the  system  are  discussed  and  portions  used  in  the  model  are 
indicated. 

The  flight  control  system  is  a  digital  fly-by-wire 
electronic  control  augmentation  system  which  uses  a  four- 
channel  parallel  network  of  computers,  electronic  circuit 
elements,  and  associated  wiring.  The  system  is  entirely 
electronic  from  pilot  controls/f eedback  sensors  to  the 
control  surface  actuators.  The  actuators  are  redundant 
electrohydraulic  servo  mechanisms  with  the  exception  of  the 
leading-edge  flap  system  which  is  a  hydraulically-powered 
rotary  mechanical  system.  Backup  mechanical  control  of  the 
stabilators  is  available,  and  the  stabilator,  aileron  and 
rudder  surfaces  have  a  backup  analog  Direct  Electric  Link 
in  the  event  of  a  total  digital  computer  failure.  Control 
surfaces  are  stabilators,  ailerons,  dual  rudders,  leading- 
ed''-'  and  trailing-edge  flaps.  The  stabilators,  leading-edge 
and  trailing-edge  flaps  are  capable  of  differential  movement. 

Pilot  inputs  are  through  a  conventional  control  stick 
and  rudder  pedal  arrangement.  Closed-loop  stability 
augmentation  is  provided  by  feedback  of  pitch,  roll,  and 
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yaw  rates,  normal  3nd  lateral  acceleration,  and  angle-of- 
attack  data.  These  feedback  signals  are  gain  scheduled  by 
air  data  and  angle  -  of - attack  information  to  tailor  the 
feedback  signals  to  the  current  flight  condition.  Cross¬ 
axis  control  signal  interconnects  are  provided  which 
improve  control  and  feedback  coordination  and  reduce  inertia 
coupling.  Reference  5  should  be  consulted  for  a  more 
in-depth  discussion  of  the  feedback  and  gain  schedule 
design  theory. 

Initially,  the  computer  model  is  to  be  used  for  studies 
involving  only  the  "up  and  away"  flight  mode,  meaning  that 
the  simulated  aircraft  is  in  a  trimmed,  stable  condition  in 
normal  flight  prior  to  the  initiation  of  a  maneuver.  This 
allows  the  model  to  be  reduced  from  the  full  control  system. 
The  assumptions  which  were  made  in  reducing  the  model  are 
discussed  in  Reference  1,  however,  several  will  be  repeated 
here : 

1.  The  aircraft  is  in  the  Auto  Flaps  Up  configuration.  The 
control  law  gain  schedules  used  are  designed  for  cruise  and 
combat  maneuvering  flight.  The  leading  and  trailing-edge 
flaps  are  automatically  positioned. 

2.  Inner  loop  control  is  being  used.  This  means  that  the 
pilot  is  the  source  of  control  inputs  through  the  stick  and 
rudders.  Outer  loop  (autopilot)  control  is  not  modeled. 

3.  The  Control  Augmentation  System  is  in  use.  This  is  the 
optimum  situation  indicating  that  the  digital  computers, 


feedback  sensors,  and  air  data,  acceleration,  and  angle-of- 
attack  measurements  are  all  operating  normally. 

4.  Trim,  external  stores,  speedbrakes,  and  anti- spin 
features  are  not  modeled.  Active  Oscillation  Control  is 
also  left  out,  as  it  is  a  function  of  external  store 
loading . 

The  model  thus  used  is  acceptable  for  simulating  the 
full  range  of  combat  maneuvering  flight  of  which  the  F/A-18 
is  capable.  Take-off  and  approach/ landing  phases  are  not 
s imulated . 

In  the  McDonnell  Aircraft  Company's  report,  F/A-18 
Flight  Control  Electronic  Set  Control  Laws  [Ref.  6],  the 
software  design  and  nomenclature  of  the  F/A-18  digital 
control  system  is  described.  Those  portions  which  are 
included  in  the  model  are  briefly  reviewed  below,  including 
for  completeness  the  longitudinal  portion  which  was 
discussed  in  Reference  1.  Readers  desiring  in-depth  infor¬ 
mation  of  the  F/A-18  flight  dynamics  and  control  system 
theory  of  operation  should  refer  to  References  S  and  6. 
LONGITUDINAL  AXIS 

The  control  system  consists  of  the  following  paths: 
o  Stick  pitch  input 
o  Pitch  rate  feedback 
o  Normal  accelerometer  feedback 
o  Angle-of-attack  feedback 
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o  Inertial  decoupling  feedback  (roll  rate  *  yaw  rate) 
o  Forward  loop  integrator 

These  paths  are  sunned  tc  form  the  stabilator  command. 

Except  for  inertial  decoupling  feedback,  all  of  the  paths 
are  used  as  inputs  to  the  forward  loop  integrator  path, 
which  reduces  the  steady-state  difference  between  maneuver 
command  and  aircraft  response  to  zero.  The  signal  paths 
are  gain  scheduled  by  air  data,  angle- of - attack ,  and 
acceleration  functions  which  are  discussed  in  Chapter  III. 

A  digital  notch  filter  is  located  after  the  feedback  summing 
junction  to  attenuate  structural  vibrations.  In  the  Auto 
Flap  Up  mode*  gain- scheduled  angle- of - attack  information  is 
used  to  provide  maneuvering  flap  commands. 

LATERAL  AXIS 

The  control  system  consists  of  the  following  paths: 
o  Stick  roll  input 
o  Roll  rate  feedback 
o  Rudder  pedal  interconnect 

The  lateral  command  is  formed  by  the  sum  of  the  signal  paths, 
which  are  gain- scheduled  similarly  to  the  longitudinal  axis 
paths.  The  stick  and  roll  rate  paths  incorporate  digital 
notch  filters  to  attenuate  structural  vibration  inputs. 
Differential  stabilator,  and  leading-edge  and  trailing  edge 

*  flap  commands  are  separately  gain- scheduled  and  directed  to 

* 

f  their  respective  paths  in  the  longitudinal  axis  control 

system. 
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DIRECTIONAL  AXIS 

The  control  system  consists  of  the  following  paths: 
o  Rudder  pedal  yaw  input 
o  Yaw  rate  feedback 
o  Lateral  acceleration  feedback 
o  Rolling  surface  interconnect 

o  Inertial  decoupling  feedback  (pitch  rate  *  roll  rate) 
The  rudder  command  is  the  sum  of  s tructural - filtered ,  gain 
scheduled  feedback  signals,  and  gain-scheduled  rudder  pedal 
and  rolling  surface  interconnect  signals. 

Actuator  data  from  Reference  5  was  not  expressed  in 
terms  of  damping  and  natural  frequency,  therefore,  data 
that  was  available  on  an  actuator  very  much  like  the  F/A-18 
stabilator  actuator  was  used  as  a  guide.  It  was  decided  to 
use  a  second  order  model  for  all  of  the  control  surface 
actuators,  and  to  use  the  same  damping  ratio  and  natural 
frequency,  0.7  and  40Hz,  respectively,  until  the  data 
specified  in  that  manner  was  obtained  for  the  other 
actuators . 

Digital  filters  other  than  the  structural  filters 
already  mentioned  include  lag,  lead-lag,  and  integrators. 
The  constants  for  all  of  the  filters  are  listed  in 
Chapter  16  of  Reference  5.  Aliasing  filters  were  treated 
as  analog  filters  (s-domain)  as  they  always  occur  prior  to 
A/D  converters.  Digital  filters  used  the  z-domain 
representation . 


14 


In  Chapter  III,  the  modeling  methodology  by  which  th 
F/A-18  digital  control  system  as  described  briefly  above 
was  coded  for  computer  study  using  CSMP  is  discussed. 


III.  PROGRAM  METHODOLOGY 


The  primary  source  of  information  used  in  making  the 
computer  model  was  the  McDonnell  Aircraft  Company  F/A-18 
Flight  Control  System  Design  Report  [Ref.  5].  Figures 
16.1,  16.2,  and  16.3  of  Reference  5  are  the  block  diagrams 
of  the  longitudinal,  lateral,  and  directional  control 
systems,  respectively.  Chapter  16  of  Reference  5  contains 
descriptions  of  system  operation,  control  law  algorithms, 
and  digital  filter  specifications.  Since  there  have  been 
several  versions  of  the  control  laws  to  date,  it  should 
be  noted  that  the  version  used  in  this  program  is  OPY  8.2.1, 
current  as  of  31  August  1982. 

Program  modeling  began  by  reducing  the  control  system 
block  diagrams  (Figures  16.1,  16.2,  16.3  of  Reference  5) 
to  the  forms  desired  for  research.  This  was  done  by 
applying  the  assumptions  listed  in  Chapter  II  and  in 
Reference  1.  Next,  scheme  of  labeling  the  control  system 
paths  was  developed.  Readers  of  Reference  1  should  note 
that  the  labeling  system  used  there  has  been  changed  in 
this  program  to  accommodate  the  more  general  nature  of  the 
full  three-axis  model.  The  paths  which  are  modeled  are: 


1. 

Pilot  inputs:  pitch, 

roll , 

and  yaw- 

2. 

Angle-of-attack  feedback 

3. 

Rate  gyro  feedbacks: 

pitch , 

roll,  and  yaw 
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4. 

Normal  accelerometer 

feedback 

5 . 

Lateral  accelerometer  feedback 

6. 

Stabilators,  symmetrical  and  differential 

/  • 

Ai lerons 

8. 

Rudders 

9. 

Leading-edge  flaps, 

symmetrical 

and  differential 

10. 

Trailing-edge  flaps, 

symmetrical 

and  differential 

The  labels  and  nomenclature  for  each  path  are  listed  below 
PILOT  INPUTS 

PP  pilot  pitch  path 

PR  pilot  roll  path 

PY  pilot  yaw  path 

ANGLE -OF- ATTACK  FEEDBACK 
AA  angle-of-attack  path 

ALPHAT  computed  angle-of-attack 
RATE  GYRO  FEEDBACKS 
PG  pitch  rate  gyro  path 

RG  roll  rate  gyro  path 

YG  yaw  rate  gyro  path 

CSAOAT  approx,  cos (ALPHAT) 

SNA OAT  approx,  sin (ALPHAT) 

P2A,B,C  filter  P2  arguments 
Y3A,B,C  filter  Y3  arguments 
NORMAL  ACCELEROMETER  FEEDBACK 
NZ  normal  accelerometer  path 

NZA  nz  for  gain  schedule  use 
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NZAF 


filtered  incremental  nz 


NZAR1  nca-based  roll  rate  gyro  parameter 

P5A,B>C  filter  P5  arguments 

LATERAL  ACCELEROMETER  FEEDBACK 

NY  lateral  accelerometer  path 

STABILATOR  PATH 

ST  main  stabilator  path 

DS  differential  stabilator  path 

ST2R  main  to  differential  stabilator  path  signal 

RST  right  stabilator  path 

RSTDEF  right  stabilator  deflection 

LST  left  stabilator  path 

LSTDEF  left  stabilator  deflection 

AILERON  PATH 

AL  main  aileron  path 

YV1R  rudder  to  roll  crossfeed  path 

RAL  right  aileron  path 

RALDEF  right  aileron  deflection 

LAL  left  aileron  path 

LALDEF  left  aileron  deflection 

RUDDER  PATH 

RD  main  rudder  path 

RSR  roll  surface  to  rudder  interconnect  path 

LRD  left  rudder  path 

LRDDEF  left  rudder  deflection 

Y5A, B,C  filter  Y5  arguments 
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LEADING -EDGE  FLAP  PATH 
LE  main  leading-edge  flap  path 

DLE  differential  leading-edge  flap  path 

RLE  right  leading-edge  flap  path 

RLEDEF  right  leading-edge  flap  deflection 

LLE  left  leading-edge  flap  path 

LLEDEF  left  leading-edge  flap  deflection 

TRAILING- EDGE  FLAP  PATH 
TE  main  trailing-edge  flap  path 

DTE  differential  trailing  edge-flap  path 

RTE  right  trailing-edge  flap  path 

RTEDEF  right  trailing-edge  flap  deflection 

LTE  left  trailing-edge  flap  path 

LTEDEF  left  trailing-edge  flap  deflection 

Nomenclature  from  the  McDonnell  Aircraft  Company  literature 
[Refs.  5  and  6]  has  been  included  where  it  was  deemed 
helpful  for  reference.  This  primarily  includes  such  terms 
as  PK_,  RK_,  YK_,  PV_,  RV_,  and  YV_;  these  are  constants  or 
variables  which  occur  at  significant  points.  The  terms 
QC  (dynamic  pressure)  and  PS  (static  pressure)  are  also 
from  the  McDonnell  Douglas  Aircraft  Company  literature. 
Carter  [Ref.  1]  lists  and  defines  some  additional  terms 
which  have  been  used  in  this  program,  including  filter  and 
constant  nomenclature.  The  remainder  of  the  terminology 
was  developed  specifically  for  this  program.  Block  diagrams 
incorporating  the  simplying  assumptions  mentioned 
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previously  and  using  the  above  notation  are  included  in 
Appendix  A.  These  particular  block  diagrams  are  intended 
only  to  aid  in  identifying  the  various  terms  with  their 
location  in  the  diagrams  and  leave  out  much  of  the  detail 
found  in  Figures  16.1,  .2,  .3  of  Reference  5. 

Carter  [Ref.  1,  Chapter  III]  gives  an  excellent  des¬ 
cription  of  the  methodology  of  using  CSMP  to  model  the 
block  diagram  elements  of  Figures  16.1,  16.2,  and  16.3 
of  Reference  5.  Some  changes  have  been  made  to  those 
procedures,  however: 

1.  The  control  law  gain  schedules,  referred  to 

singularly  as  Function  F _ ,  have  been  removed  from  nosort 

sections.  Instead,  all  of  the  gain  schedules  have  been 
coded  as  Fortran  functions  and  called  as  required.  This 
step  has  removed  over  300  lines  of  code  from  the  CSMP 
program,  since  Fortran  subroutines  or  functions  do  not 
count  against  the  allowable  number  of  CSMP  statements. 
Twenty- four  functions  have  been  added  for  the  lateral  and 
directional  axes,  in  addition  to  13  for  the  longitudinal 
axis  discussed  in  Reference  1: 

FUNCTION  DESCRIPTION 

F4  Roll  Rate  Feedback  Gain  Schedule  (QC,PS) 

F6  Differential  Stabilator  Gain  Schedule  (RI,PS,RV11) 

F7  Lateral  Command  Gain  Schedule  (QC,PS) 

F10  Rudder  Command  Gain  (QC) 

F13  Lateral  Command  Gain  Schedule  ( RI , PS .STORES) 

20 


I 


F17  Rudder  Pedal  Command  Gain  Increment  (AOA) 

F30  RSRI  Gain  Schedule  (QC,PS) 

F31  Differential  T.E.  Flap  Gain  Schedule  (RI,PS) 

F34  Differential  T.E.  Flap  Gain  Schedule  (AOA) 

F3S  Lateral  Forward  Loop  Gain  Schedule  (AOA) 

F36  Aileron  Gain  Schedule  (QC, PS, RI) 

F38  RSRI  Gain  Schedule  (AOA,RI,PS) 

F39  Rudder-Roll  Interconnect  Gain  Schedule  (AOA) 

F41  Roll  Surface  Limit  Schedule  (AOA,RI) 

F42  RSRI  Nonlinear  Gradient 

F45  Directional  Forward  Loop  Gain  Schedule  (QC,PS) 

F90  Lateral  Acceleration  Feedback  Gain  Schedule  (RI,PS) 

F93  Differential  L.E.  Flap  Gain  Schedule  (RI,PS,NZAF) 

F96  Yaw  Rate  Gain  Schedule  (QC,PS) 

F101  Differential  Stabilator  Load  Alleviation  (RI,PS, 

NZAF) 

F108  Directional  Inertial  Gain  Schedule  (QC) 

F112  Lateral  Acceleration  Gain  (RI) 

F113  Lateral  Acceleration  Gain  (AOA) 

F114  Rudder  Pedal  Command  Gain  Increment  (PI) 

Algorithms  for  the  implementaion  of  the  functions  are  found 
on  Chapter  16  of  Reference  5. 

2.  Frequency  averagers  and  rate  limiters  are  modeled 
using  CSMP  macro  statements.  Thus,  all  nosort  sections  have 
been  eliminated  from  the  current  3-axis  program. 
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3.  Cross-axis  signal  paths  for  differential  stabilatcr 
leading-edge  and  trailing-edge  flaps  have  been  added. 
Rolling-surface- to- rudder  interconnect  (RSRI)  ,  rudder-to- 
roll  crossfeed,  and  inter-axis  distribution  of  rate  gyro 
feedback  signals  have  been  incorporated. 

4.  As  noted  in  Reference  1,  in  the  longitudinal-only 
simulation,  the  stabilators  and  flaps  were  modeled  for  one 
direction  of  motion  only;  this  constraint  has  been  removed 
in  the  current  3-axis  program. 

Descriptions  of  the  algorithms  used  to  implement  A/D 
and  D/A  converters,  digital  filters,  aliasing  filters, 
limit -functions ,  and  actuator  servomechanisms  are  contained 
in  Reference  1.  They  are  used  in  the  current  program 
without  modification. 
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IV.  MODEL  VALIDATION 


The  process  of  validating  the  computer  model  of  the 
F/A-18  flight  control  system  was  carried  out  in  two  major 
steps  : 

1.  Individual  block  diagram  elements  were  tested  using 
probable  ranges  of  input  values.  Digital  filter  macros, 
rate  limiters,  frequency  averagers,  and  gain  schedule 
functions  are  included  here.  Carter  [Ref.  1]  describes 
the  testing  of  all  of  the  block  diagram  elements  mentioned 
with  the  exception  of  the  gain  schedule  functions  relevant 
to  the  lateral  and  directional  axes,  which  were  individually 
verified  by  this  researcher  exactly  as  were  the  gain 
schedule  functions  used  in  Reference  1. 

2.  The  entire  model  was  assembled  and  subjected  to 
inputs  from  the  stick,  rudder,  and  angle- of - attack , 
acceleration,  and  feedback  sensors.  The  goal  of  this  step 
was  to  verify  correct  direction  of  motion  of  control 
surfaces  in  response  to  unambigous  inputs. 

Inasmuch  as  Reference  1  has  already  discussed  the  procedures 
and  results  relevant  to  Step  1,  with  the  exceptions  noted, 
the  current  discussion  will  involve  Step  2.  Further, 
since  Reference  1  contains  the  validation  results  for 
stabilator,  leading-edge  and  trailing-edge  flaps  due  to 
pitch  inputs,  here  the  control  surface  motion  due  to  roll 
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and  yaw  inputs  will  be  of  primary  concern.  The  source 
of  information  used  as  a  reference  in  comparing  model 
performance  to  design  specificat ion  was  the  McDonnell 
Aircraft  Company  F/A-18  Flight  Control  System  Design 
Report  [Ref.  5] . 

The  computer  model  requires  the  following  inputs: 

1.  Static  pressure  (PS)  in  Ib/sq.ft. 

2.  Dynamic  pressure  (QC)  in  lb/sq.ft. 

3.  Pilot  pitch  input  (PP1)  in  lb.  +  *  aft  stick. 

4.  Pilot  roll  input  (PR1)  in  lb,  +  =  right  stick. 

5.  Pilot  yaw  input  (PY1)  in  lb,  +  *  right  rudder. 

6.  Pitch  rate  gyro  feedback  (PG1)  in  deg/sec, 

+  =  nose -up. 

7.  Roll  rate  gyro  feedback  (RG1)  in  deg/sec, 

+  =  right  roll. 

8.  Yaw  rate  gyro  feedback  (YG1)  in  deg/sec, 

+  =  right  yaw. 

9.  Angle-of-attack  (AA1)  in  deg. 

10.  Normal  acceleration  (NCI)  in  "g",  +  =  nose-up  motion. 

11.  Lateral  acceleration  (NY1)  in  "g" ,  +  =  nose-right 
motion . 

Depending  on  the  need,  inputs  can  be  programmed  in  CSMP  as 
step  or  ramp  functions  (see  Reference  4) ,  or  using  any 
standard  Fortran  function  such  as  SIN  or  EXP,  for  example. 
Combinations  thereof  are  acceptable,  as  well.  An  example 
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of  a  test  roll  input  consisting  of  a  61b  right  stick  input 
for  two  seconds  followed  by  a  61b  left  stick  input  as 
written  in  CSMP  would  be: 

PR1  =  6 . 0 *STEP (0.0)  -  12  .0*STEP(2  .0) 

For  this  purposes  of  the  tests  described  in  this  thesis, 

PS  and  QC  are  held  constant  through  test  maneuvers,  though 
there  exists  no  such  constraint  in  the  model.  The  effects 
of  varying  PS  and  QC  were  verified  to  be  correct  when  each 
gain  schedule  function  was  compared  to  the  graphical  data 
in  Chapter  16  of  Reference  5.  Rate  gyro  feedback,  angle- 
of-attack,  and  acceleration  inputs  may  currently  be  pro¬ 
grammed  as  described  for  pilot  inputs.  However,  accurate 
representation  of  these  latter  three  types  of  inputs  will 
not  be  possible  until  the  F/A-1S  aerodynamic  build-up  by 
Raithel  [Ref.  2]  and  the  non-linear  equations  of  motion 
are  incorporated  into  the  program.  As  that  time,  those 
inputs  will  be  determined  by  the  program  and  will  not  be 
explicity  stated.  For  the  tests  described  here,  rate  gyro 
feedback,  angle-of-attack,  and  acceleration  data,  when 
necessary  to  simulate  aircraft  responses,  will  be  provided 
by  "best  guess"  estimation  of  those  data,  given  the  pilot 
input . 

Figures  presented  in  this  chapter  are  based  on  a  time 
scale  of  four  seconds  per  maneuver.  Control  surface 
deflections  are  in  degrees.  Static  pressure  and  dynamic 
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pressure  inputs  are  listed  on  each  figure  in  terras  of 
altitude  and  Mach  number.  Standard  day  conditions  are 
assumed  in  all  cases.  Angle  -  of - attack  is  occasionally 
listed  where  it  plays  a  significant  role  in  shaping  the 
control  system  response. 

Figure  4.1.  depicts  differential  stabilator  response  to 
a  +/-  6.01b  lateral  stick  input.  The  shift  from  right  to 
left  stick  takes  place  at  the  2.0  second  mark.  The  per¬ 
tinent  gain  schedule  functions  are  Functions  6  and  101. 

For  the  given  flight  conditions  and  pilot  input,  the 
stabilator  response  is  very  close  to  maximum.  The  response 
is  decreased  at  higher  dynamic  pressures. 


'o.oo  0.«0  0.80  l.»  1.6C  2.00  2.%0  2.80  3.20  3.8C  4.00 

TIME 


Figure  4.1.  Differential  Stabilator  Response  to  Roll 


Figure  4.2.  shows  differential  stabilator  response  to 
the  same  input  as  for  Figure  4.1.,  but  with  roll  rate  gyro 
feedback  added  as  a  ramp  type  input.  This  causes  the  damp¬ 
ing  effect  on  the  stabilator  motion.  The  overshoot  caused 
when  the  input  signals  are  reversed  in  sign  is  to  some 
extent  believed  to  be  an  artifact  of  the  programming  of 
the  roll  rate  feedback  input,  which  only  approximates  an 
actual  feedback  input.  This  phenomena  can  be  noticed  on 
some  of  the  following  graphs,  as  well. 

Figure  4.3.  shows  differential  stabilator  motion  due  to 
a  rudder  input.  This  demonstrates  the  rudder- to- roll 
surface  CAS  interconnect,  which  is  gain-scheduled  by 
Function  39  using  angle-of -attack.  A  ramp  input  was  used 
to  simulate  the  angle-of-attack  signal. 

Aileron  motion  due  to  a  +/-  6.01b  roll  input  is 
depicted  in  Figure  4.4.  Aileron  gain  schedules  are 
Functions  35  and  36,  which  accept  angle-of-attack  and  air 
data  inputs,  respectively.  Both  functions  decrease  gain 
as  their  input  values  increase:  in  the  case  of  increasing 
angle-of-attack,  this  is  to  reduce  sideslip;  the  air  data 
gain  schedule  lessens  aileron  response  at  high  dynamic 
pressures  where  aeroelastic  effects  would  otherwise  cause 
aileron  reversal. 
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Figure  4.2.  Differential  Stabilator  Response  with 
Roll  Rate  Feedback 

When  roll  rate  feedback  is  added  to  the  aileron  roll 
motion, it  is  seen  that  a  damping  effect  is  present,  as  it 
should  be.  This  is  shown  in  Figure  4.5.  Roll  rate  feedback 
is  gain-scheduled  by  Function  4,  which  acts  to  decrease 
feedback  gain  as  dynamic  pressure  increases.  The  overall 
lateral  control  system  gain  is  shaped  additionally  by 
Functions  7  and  13,  using  air  data  inputs.  These  gain- 
schedules  do  not  easily  lend  themselves  to  straight-forward 
explanations  of  purpose,  indeed,  there  are  certainly  multi¬ 
ple  purposes  for  the  design  of  these  functions.  Reference  5 
is  more  explicit  as  concerns  the  gain-schedule  design 
philosophy . 
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The  rudder- to - rol 1  surface  CAS  interconnect  signa 


causes  an  aileron  deflection  as  angle-of -attack  is 
increased.  This  effect  is  shown  in  Figure  4.6.,  and  is 
analogous  to  that  noted  in  Figure  4.3. 

At  higher  dynamic  pressures,  as  noted  earlier,  aileron 
gain  is  decreased  to  reduce  aeroelastic  effects.  Differen¬ 
tial  leading-edge  and  trailing-edge  flaps  are  incorporated 
to  maintain  an  acceptable  roll  response  under  these  cir¬ 
cumstances.  Differential  leading-edge  flap  motion  in 
response  to  a  +/-  10.0  lb  roll  input  is  displayed  in 
Figure  4.7.  At  the  indicated  Mach  number,  aileron  response 
is  nil.  Function  93  governs  the  response  as  a  function  of 
air  data  and  normal  acceleration  inputs. 


Figure  4.7.  Differential  Leading-Edge  Flap  Response  to  Roll 
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Figure  4.8.  Differential  Leading-Edge  Flap  Response  to  Roll 
with  Roll  Rate  Feedback 

The  expected  effect  of  roll  rate  feedback  on  differen¬ 
tial  leading-edge  flap  motion  is  shown  in  Figure  4.8.  The 
damping  is  noticeably  less  than  for  the  differential 
stabilator  and  aileron  damping  observed  previously;  this  is 
due  to  the  inherently  higher  damping  of  aircraft  rolling 
motion  at  the  dynamic  pressures  where  differential  flap 
motion  is  required.  Figures  4.9.  and  4.10.  depict  similar 
responses  for  differential  trailing- edge  flaps,  without  and 
with  roll  rate  feedback,  respectively.  Functions  31  (air 
data)  and  34  (angle-of- attack)  are  the  gain  schedules 
responsible  for  differential  trailing-edge  flap  motion. 
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Figure  4.9.  Differential  ‘Trailing-Edge  Flap  Response  to  Roll 
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Figure  4.10.  Differential  Trailing-Edge  Flap  Response  to  Roll 
with  Roll  Rate  Feedback 
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Figure  4.11.  Rudder  Response  to  Yaw 

Directional  control  system  responses  are  made  by  the 
twin  rudders,  which  do  not  act  differentially  in  the  Auto 
Flap  Up  mode.  Rudder  response  to  a  +/-  20. 01b  rudder  pedal 
input  is  shown  in  Figure  4.11.  The  gain  schedules  which 
shape  the  motion  are  Functions  10  (air  data)  ,  71  (air  data)  , 
and  114  (angle-of-attack) .  Figure  4.12.  shows  rudder  motion 
when  feedback  signals  are  added  to  the  simulation.  The 
damping  effect  is  the  sum  of  stability  axis  yaw  rate  (this 
feedback  signal  is  a  blend  of  yaw  and  roll  rates,  and  angle- 
of-attack  data),  and  lateral  acceleration  feedbacks. 
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Figure  4.12. 


Rudder  Response  to  Yaw  with  Feedback 
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Figure  4.13. 


Rudder  Response  to  Roll 
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Several  gain  schedules  are  involved  here  in  shaping 
rudder  response  to  the  requirements  of  sideslip  reduction, 
prevention  of  excessive  vertical  tail  loads,  inertia 
coupling  reduction,  and  high  angle  -  of- attack  maneuvering 
stability.  Reference  5  discusses  these  considerations. 

Rolling  surface-to- rudde r  interconnect  signal  operation 
can  be  seen  in  Figure  4.13.  Functions  30  (air  data),  and 
38  (angle-of -attack  and  air  data)  shape  this  response. 

Verification  of  all  of  the  signal  paths  of  the  F/A-1S 
flight  control  system  in  the  Auto  Flap  Up  mode  has  thus 
been  accomplished.  The  gain  schedules  (Functions)  were 
each  individually  verified  prior  to  incorporation  into  the 
computer  model . 

The  process  of  model  validation  has  thus  been  to 
observe  model  responses  to  some  relatively  basic  inputs, 
and  then  to  judge  these  responses  as  being  plausible  or  not 
with  respect  to  the  data  provided  in  Reference  5.  The 
computer  model  has  given  logical  results  as  noted  in  this 
chapter,  and  also  for  a  wider  range  of  input  signals  and 
flight  conditions  than  have  been  discussed  here.  The 
conclusion  is  that  the  flight  control  system  model  is  valid 
representation  of  the  F/A-18  flight  control  system  within 
the  range  of  the  simplifying  assumptions  and  limitations 
noted  previously  in  Reference  1  and  in  this  thesis. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  computer  simulation  of  the  F/A-18  digital  electronic 
flight  control  system  developed  here  is  suitable  for  further 
use  in  aircraft  control  systems  studies.  With  the  inclusion 
of  F/A-18  aerodynamics  and  equations  of  motion,  it  will  also 
serve  to  examine  simulated  maneuvering  flight  under  condi¬ 
tions  found  near  or  at  the  accepted  flight  envelope.  Com¬ 
parative  studies  may  be  done  of  newer  control  systems 
concepts,  using  the  known  response  as  the  basis  for 
comparison . 

CSMP  has  size  limitations  for  several  internal  para¬ 
meters,  the  most  pertinent  of  which  are  total  number  of 
program  statements  (1900)  ,  and  maximum  number  of  statements 
in  any  given  sort  section  (600)  .  The  current  program  is 
right  at  the  latter  limit  with  599  statements  in  one  sort 
section.  Some  statements  could  be  combined  to  reduce  this 
number;  this  is  not  recommended  until  the  user  is  familiar 
with  the  F/A-18  flight  control  system  and  its  representation 
here.  It  is  recommended  that  the  aerodynamics  and  equations 
of  motion  be  written  as  Fortran  subroutines  and  called  at 
the  beginning  of  the  dynamic  section.  This  would  use  only 
two  lines,  thus  only  two  additional  statements  would  have 
to  be  combined  to  make  the  necessary  room.  There  are  not 


obvious  dividing  lines  where  the  current  program  could  be 
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broken  up  into  more  than  one  sort  section,  however,  users 
experienced  in  CSMP  may  find  this  to  be  possible  when  the 
aerodynamics  and  equations  of  motion  are  incorporated. 
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APPENDIX  A:  MODEL  BLOCK  DIAGRAMS 


PILOT  PITCH  INPUT  PATH 
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ANGLE-OF-ATTACK  SENSOR  PATH 


STA8ILAT0R  PATH  (CONT.) 
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AILERON  PATH  (CONT.) 


LEAO ING-EDGE  FLAP  PATH  (CONT.) 
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APPENDIX  3: 


MODEL  COMPUTER  PROGRAM 


******  ****************************  ***************** 

**  f-18  FLIGHT  CONTROL  LAN  SISULATION  ** 

**  WW 
****  **********  ********************  ***************** 

********************  *******  ************************ 

*****INTEGRATISG,IAG,AND  NOTCH  DIGITAL  FILTERS*.*.* 

4  m**  **************  ***  *****  *******  ***********  **  *  ** 

SAC  EC  FC0T*ZINT  (FIN  ,  K 1 , 8E  .  IBP ,  POOTZ 1) 

PROCEDURAL  „  „„ 

IF  (IBP.  N E.  1.0)  GO  TO  10 
IF  (KEEP. NE. 1.0)  GO  10  10 
IF  (TIME.  EC»  0.  0)  GO  TO  10 
FO  U1*KA*  FI N  ♦  KB*FCCTZ 1 
FOOTZ1-FOUT 

10  CONTINUE 

S  AcloSoDT-ZLAG  (FIN  ,  HA ,  KB  ,  KC,I  HP#FI  XZl  ,  FOD  TZ 1 ) 

PROCEDURAL 

IF  (ISP.N E.  1.01  GC  1C  JO 
IF  (KEEP.  NE.  1.0)  GC  TO  10 
IFjTISE.  EO.C.O)  GC  10  10 
FOUT-KA*  FIN  -  KE*FINZ 1  ♦  KC*FOUTZ1 
FO  UTZ1 *FO  U1 
FI  NZ  1=FI  N 

10  CONTINUE 

IacROCFOOT=ZN01CH(FIS,KA#KE,KC/KD,KE,IHE,FINZ1,FINZ2.FCUTZ1,FODTZ2) 

PROCEDURAL 


10 


ENDBAC 


IF  (ISP.NE.  1.0!  GC  10  10 
IF  (KEEP.  NE.  1.  0)  GC  10  10 
IF  IlIHE.  EQ  .C.O)  GO  1C  10 
FOUT=KA*FlN  +  kE*FINZ1 
FOUTZ2  =  FCCTZ1 
FOOTZ1  *  FCOT 
FINZ2  *  FINZ1 
FI NZ 1  *  FIS 
CON1INUE 


+  KC*FINZ2  -  KI) *FO DTZ 1  -  K3*FODTZ2 


****** ***•******••••*****»••* 

•  *»**F6EC0ENCI  IV  F SAGERS  ***** 
***************************** 


TO  20 
TO  5 


MACRO  Z40*A V2C«Q(Z2C,IEF) 
IHCCECDRAL  _  ,  „ 

IF  (KEEP.  NE.  1.0  GO 
IF  (TIME.NE.O.C)  GO 
Z20Z1=Z20 
Z4CZ1= Z2C 
CEI=0. 0 
GO  TO  10 

IF  (ISP.  EC.  1.0)  GC  TO  10 
Z«C*Z“ OZ  1+DSL 
GC  TO  15 
Z40*Z20 


10 


DEt=  (Z  20-Z20Z  1J/2.0 
Z2CZ1=Z2C 


ZU0Z1*Z40 

CONTINUE 

CONTINUE 


GO  TO  20 
GC  TO  5 


15 

20 

END MAC 

SACSC  Z80-AV4080  (Z40,  ISP) 
PROCEDURAL 

IF  (KEEP.NE.  1.0) 

IF  (TIME.NE.O.C) 

Z4QZ 1*  Z4  0 

zeczi*z4C 

DEI*0. 0 
GO  TO  10 

5  IF  (IMP.  EC.  1.0)  GC  TC  10 

ZEC*Z80Z1*DEL 
GC  TO  15 
10  Z€0*Z40 
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:§{ 


GO  TO 
GO  TO 


GC  TC  10 


DSL*(Z4C-Z40Z1)/2.0 
Z40Z1*  Z4  0 

zeczi-zac 

15  CCSTIHOE 

20  CONTINUE 
E  ND  B  AC 

MACHO  Z8C-AV20eC(Z2C.ISE) 

IfiCCEOUSAL 

IP  (KEEP.SE.  1 
IP  Itibe.se.  0 
Z2CZ1-Z2C 
Z80Z1-Z20 
GO  TO  10 

5  IP  (ISP.  PC.  1.0) 

Z60-Z80Z1+DEL2 
GC  TO  IS 

io  zeo»z2 o 

DEL2=(Z20-Z20Z1)/4.0 
15  CONTIS OE 

Z20Z1-Z20 
Z8CZ1-Z8C 
20  CCNTISOE 

END8AC 

***************************** 

*******  HATE  limits?  ******** 
****** **»•***»»*»** ********** 

8  AC  EC  SCOT- RAT  LB  (HIS, IBP,  INC) 

PROCEDURAL 

IF  (KEEP.  SE.  1)  GC  TO  100 
IF  (USE. SE. 0.0)  GO  TO  80 
60UTZ1-BIN 
80  CONTINUE 

CEL  1*  SIN -SCOTS  1 
0ELII3*LI MIT (-ISC. I  SC, DELI ) 
SCUTA-BOUTZI+DELLlS 
IF  (IBP. SE.  1.0)  GC  IC  90 
SOUTZ1-BOOTA 
90  CONTINUE 

SOOT-ZHOID  (IMP.  SCOT  1) 

100  CONTINUE 

ENDBAC 

• 

INITIAL 

* 

*  **•*••**«*•*•*•**********•**••*•* 

*  DIGITAL  FILTER  INITIAL  CONDITIONS 

*  ***•••••**•*•*•******•*•***•*•*** 


PS.QC.PK.BK.YK  ASC  HATE  CONST  APTS 


CONSTANT  PS-1761. O.CC-800.0, PK9— 1  .  1543, P K10-0. 46  47 , PK12-0. 5654, . 
BK  1-  1.  177  .RK2-3.  ”.FK7-50.Q  ,  BK  10-1.  33,  HK  11-  190.0,  3K13-1.0,..  . 
IK7-.S1271.IK8-.  046  77.  TK1S- 1.31266.  rK16  -.024  69,  ..  . 

BATE20-0.  05  ,  SAT  E4C-0 . 3  25,  BATES  0-0.0  125 


DYNAMIC 

* 

SCST 


♦••••INPUTS***** 


PP 1-00 . 0*STEP  (C.O) 
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CJO  II 


**  *  * 


PE  1=0.  C*  STEP  (C.C)  -00  .0*STEP(2.C) 

EY 1 =20 . 0  *  STEP  (C.0) -40.0*  STEP  (2.0) 

EG1=00.0*S?£?  0.0) 

RG1=00.0»EA«P  <C.  C  -00.0*RAHP  (2.0) 

YG  1=0.  5 * E IMP  (0.0)  -1.  5*HAHP  (2. 0) 

AA1=7.  0* STEP  (0.0  ♦  0.C*SASE{0.  CJ 
SZ1=0. C*SIEP (0.0) 

NY  1=0.  0*  EASE  (0.0)  -0.  0*E  AMP  (  2  .  Cl 

BI=QC/PS 

FIC*tI»IT  (0.  0,  1800.0, PS) 

QKF=LI3IT  (200.C,2C00.0,«JC) 

*************************** 

*****r.YP0IS2  FCNC7IONS***** 

*************************** 

IKE20=  I3E0LS  (0.0.EATI20) 

IMP40* iaEOLS (0.0.SATE40 
IBP80=iaECLS (C.0,EAIE80) 

**•••• ******************************* 

*****A3GLE  OF  ATTACK  SEHSCB  PATH***** 

******  *********************  ********** 

AJ2=C3PXFI/0.  0,0.0,0.74,2  09. 0,AA  1*43  63  1.0) 

AA3*2HOLG  (XMPaQ,  AA2) 

ALPHAT=0.59*AA3*1 .9 
A  A6=  AX  FH  AT -22.0 
AA7=LTMIT  (0. 0,1CCCO.C.AA61 
EK17=ZHCld(XHE20,E37  (32 AJ  *F3  2A(QKF)  ) 

P»1=PK 1 7*AA7 
E71A*ZHC1E(I3E20,E71) 

******  *******  *******  ***** 

*****3ATE  GYRC  PA IH  S  ***** 

I************************ 

EG2=C3?XE1(0.  0,0.0,0.89,7  8.5,  EG1*6 16  2. 25) 

PG3=ZHOLC  <I3P80,FG2) 

PK11  =  1.65*F22  ICC) 

E2A=1. 0  +  EK11*  (1.U-EK121 
P2E»<1  .0+EK1 1)  *  (1.0-EK12) 

E2C=  1.  0-EK12 

PGU=ZLAGJEG3,E2A,E2E,P2C,ISB80,PG3Z1 ,PG4Z1) 

E7U=I.IHIT  (-80.0. 120.  fl,?G4) 

PGS=PV4*  (  (F6  8  (CC)  *F3  2A  (QKF)  )  +fUQ  (HI ,  E5 ,  CC,  PIQ)  ) 
F74A=ZHOIC(I3E20' EV4) 

RG2=C.lPXEljO.C,Q.C,0.8,90.0,3G1*  8100.0) 

RG3=ZHOLE (iaPSC, SG2) 

BG4=z;i0TCH(BG3,.  16  5  1 8, .  3303  6 ,  .  16  5 1 8,  - .  5  0084,  .  1 6 1 57 ,  IHP80  ,  .  .  . 
SG3Z1.SG3Z2,3G4Z1 ,3GUZ2) 

8G5=ZH  OT  CBJRG4,.  7589  2,-  1.  24§  3  1 , .  7  D53  4 ,-  1 . 248  31  ,  .  4  64  26 , 13P8  0  , .  . 
3G4FZ1  ,EG4FZ2  ,  EG5Z1  ,FGSZ2) 

BG6=ZMOTCH(HG5,.6  1223,-.  275  05  ,  .  5  775  1  ,- .  27505,  .  1  8974  ,  IHP80 ,  ... 

RG5?Z1,PG5?Z2,RG6Z1 ,EG6Z2) 

BG6A=DEAXSP(-RVlS,BY  16,  SG6) 

BK6=F4  (QC.PS) 

BG6B=LISIT(-10.0,  10 . 0, 3G6)  *LI8IT  (0.0, 0.  12,  .12-3K6) 


BG6  B=LISIT  {  - 10 . 0 ,  10 . 0,BG6)  «LISIT  (0.  0 , 0.  12,  .12-3K6) 
BG7*RG6*BK6*3G6A*FG6E 

YG2=C3?XEI(0. 0,0.0,0.89,  200.0,40  000.  0*YG1) 
YG3*ZHOiE(iap40,TG*) 

C£AOAT*1 .C-0. 5»< (Ai?HAT*0. 0174533)** 2.0) 
IG4=YG3*CSAOAT 


SNAOAT* (1.0-0. 125* ( (ALPHAT’O. 0174533  ) **2.0))  *AL?HAT*0. 0174 533 
IGS=YG4-SHA0AT*ZHCIC  jIHP40f  RG2) 

IG6=HATLS  <YG5,I8PU0.  15.  0 » HAT E 4(3) 

Y3A«<YK15*(1 .5-1 .S87S*YK15) *YR16) 

Y3E*YK  15*  (1.0-YK16) 


Y3C-1.0-YK16 

YG7*ZIAG  (YG6,  T3i,Y3E,Y3C,  X3P4C,YG6Z  1  ,TG7Z1) 
IG8«RAILB  (YG7  ,I8P40, 15.0*SAIE40) 


*****NCR8AL  ACCELEB08STER  PATH***** 
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#*#•**  ***** 


KZ2=C3PXEt(0. 0 .0.0,0.89, 2  00.  O.SZ  1*4  0  000.0) 
NZ3=ZH0LE  (IMPU5,SZ2) 

S  ZA  =NZ  3 

NZJH1  =  -50.0*LIKIT  (0 . 0 ,4 . 0 . SZ A)  +280.0 
3  V  16=A  MI N 1  (NZ&51.5K1  if 
SZ9=)1Z3-  (3  j3*»  2)  *  (6.SS29E-61 
NZ4A=ZH0LD (IMP  20 ,  SZ4) 

37AF=ZIST  (NZ4  A.0.25, 0.7  5  ,  IMP2  0.SZAFZ  1) 

PS  4=  (1. 0+SK9*  (1.0-EK  10)  ) 

P  53  =  (1  . 0+EK9)  *  (1.C-EK10) 

E5C= i . 0-EK10 

HZS=ZLA3  (SZ4  .E5A  .ESE  ,?5C  ,IBP4  0,SZ4Z1  ,HZ5Z1) 
EV2=  LIMIT  (-10.  C,  10.  0,:iZ5) 

3Z7=PV2*2.5*F32A (CKF) 

EV2A=zaoLt(iME20,mj 

******  *******  ********** **•*•***•*• 

****LAT£BAL  ACCiLxBQflETZS  PATH**** 

******  ****** ***************  ******* 


NX2=CMPXEL(Q.C.C.C, 0.89, 7 8.  54,6168.  53*3X1) 

sx3=zaotc  (ia ?a <5. 3X2) 

NX4=3Y  3+  (ZHOLE  <I3E40,3G3}  *(  (  1.0/347.  0)  )  »*  3 
SX5=P9C (BI,?S)*LI2IT  (-1 .  8.1.8, ST4) 

SX6  =  2I.vr  («f5,0.4,0.6,l.'!?40<  MYfiZI) 
3X7=3XS+NY6*P  1 12  (FI)  *F  1  1 3  (ALP HAT) 

••*••********»****•*»•*»* »*•**•• 

********  pilot  ia  pot  path******** 


2.3)  *  ZHOLE  (IMP40,  RG3) 


EE2* 106 .47) 


EE6=??5*  (7.0+  (0.  2* AES  (EPS))  ) 

PP7=LtMII  (-2  5.C.SC.  C.FP6) 

EV3=??7 

EV3A=ZHOLC(laE20,EV3) 

PF2=OS  AOSF (-2.C»2.0,EH1) 

EB3=CMPXEI  (0.  C.C.0.0 .7,  24  .3 ,  P32*  16 1 . 04) 

PH4=LI  M1 1  (-3.C.3.C.PF.3) 

FFS’BI  AL EL  (0.  0.0. C079,PH4) 

PR6  =  ZHCLE  (IBP80,EES) 

PB7=0ZAESE  (-0.C2  5  ,0.  C25,  PE6) 

EF3  =  P37*  ( ASS  (FF7)  *FX  1+RK2) 

?S9=ZI1QTCH  (P  36,. 6555  , -1.5512, .8234, -1.5512,. 6789, IMP80,. 

PH8Z  1,  PF8Z2.PS  9Z1,PBSZ2) 
PI2*9BADSP(-7.C.7.0.Py1  ) 

ET3*CM  ?XEI(0. 0,0. 0,1.5,  66  .6 , 4  4  25 . 6*P  12) 
?T4=ZHOLE(IMF4C,F13) 

FT5*DEADSE (-1.0,1. 0,?T4) 

YK9=0.  00072*F  10  ( QC) 

TK  10=0.234**1  C  (CC) 

FT6=?Y5*  (YK9*AES  (ET5)  +TK101 
Xm=?Y6 

EX7=  PY  6*  (0.  5-  (E  17  (AIPHATl  *F114(BI))) 

P  Y8=LI  SI  1  (-3Q.C,3C.0,?Y7) 


*****STAEILATC3  PATH***** 


FK19*ZHOl£(IMF20,J12  (31,  PS)  *?22A  (QKF1  *3.5) 
PK21*ZH01B  (I BE 20,  J  12  (HI,  PS)  *F3  2A  (QKF1  *F68  (QC)  ) 
PX22-ZHOIE [I3E20,?12  (BI  ,PS)  *F3  2A  (QKF  ) 

ST1=P7  3A*EK22 
SI2*(P74A*?K21)-ST1 
ST3=(PK19*9V  2  Al  +  P7  1 A+ST  2 
ST4-LISIT  (-10CCC.C.0.0.ST3) 

STS*LIMIT  (0.  0  ,  10000.0, aT3) 

ST6*(ST4*F23  (UPKAT)  1+ST5 
ST9-IMPL  (0.  0 , 0 . 05 ,ST 7) 

ST6A*SI6-C.5*ST8 

ST7*Zi:iT  (ST6A.0.05,  1.0,  IMP20 ,  ST7Z1) 


ST9*AV  20 4C  (LIMIT  (-50.0,  25  .0,STS)  ,IKP  20) 

ST  10=  H  Z7  +ST9 
S111=5T1C+PV1 

ST12A=ST  11+ (  (BG2*YG3) *F 107<QC)  l 
ST12=ZHOtD(THE40,£T12A) 

ST13*AV4060(ST12,ISE40) 

ST14=ST1 3+P95 

ST15=LIMI1(-25.0,25.C,ST14) 

ST  16  =  Z  SOTCH (ST15, .690 84 ,-.99066, .663  1 2, -. 99068 , . 35 3 96 , . 

IMP6C,ST15Z1.ST15Z2,ST16Z1,ST16Z2) 

ST 17=S T 1 6- (PV 3 » F 2  * A  (CXF)  ) 

ST18=Liail]-2  9.0 . 10.5.ST17) 

DS  1=  EV  2+  BV 1  A*  1  .33 
DS2=LIMIT  (-50 . 0 , 50.  0  ,  OS  1) 

DS3=DS2*EV7 

ST231*  17. 25- AES  (ST2  +  6.75) 

ST2E2=LIMIT(-5C.3.F41  (ALPHA?,  SI)  *10.  0.ST2B1) 

CS4=LI MIT  ( -aT  2 E  2 ,  ;T 2E 2,  DS  3) 

DS5=RA  TIM  (DS4,  IMP  60, 25.0*  BATE  60) 

B£T1=5T1 6-DS5 

BST2=LI MIT (-24.0,  10.  5.EST  1) 

B  ST3=Q NT2E  (8  A 1 E80  ,BS  12) 

HSSDEF=CMEXFL  (0 . 0 , 0 . 0, 0 . 7,4  0 .  0.BST3*  16  0C.0) 

L  £1 1=S 11 6+DS5 


IST3*  16CO.O) 


LST2=LIMIT(-24.0. 10.  5.LST1) 
LS13=Q!<TZ3(BATF30,LST2)i 
LSTD2?  =  CME'(PL  (0.0,0.0,5.7,40.0, 

***••*•***•******•»• 

****AILFSCN  PATH**** 

******************** 

EV 2A=  P  39  *  (  (F4  <CC*  ESI  +F13  (8I,PS))  *F7  (QC,  ES>  )  -8G7 
8V2=LIMIT  (-50.a.‘c.0.RV2A) 
avi1=F31  (SI,  PS)  lF 24  (ALPHA!) 

BK17=P6(EI,?S,5V11) 

SV7=A.1IN  1  (3X1  7 *?3 5  <AI?S AT),  BK17-F101  |SI,PS,SZA?)  ) 
TV  1F.=  LIMIT  (0.  0,1. 0,5  29  (ALPS  AT)  )*YV1 
3V1=LISII  (-200.0,200.0.  t  V  18) 
aV1A=AV4C£0(3V1,:SE40) 

All*  5V2+B71A*  1.33 
AL2*LISI1  (-50 .C,5C.0,AL1) 

AL3=AL2*E35(ALEfiAT)  **36  (Bl,QC,PS) 

AL4=LIMIT (-50.C*F41  (ALPHA!, SI)  ,50.0* 

AI5=  DEADSE (-4  2.Q,4  2.0,AL4)  +  ALU 
E All* LIMIT  (-4 2.0,25. 0,  AL5) 

8AL2=8  ATIM(F.AL1  ,IIE8C.1  00  .0*EATE80) 


F41  (ALPHA!, BI) ,AL3)*0.5 


,v .  7 ,40 .  C,-  1600 . 0*BAL3) 

A.  ill*  1-biaAA  l-4w  •  **  i.  m  0  A  L5  j 

LAI2=BATia|LAI1,IRP80, 1 00 .0*8 AIE80) 

LAL3=QNT2E(BATE80,LAL2) 

I AI02?=CSFX?L (0.0,0. 0,0. 7,40. 0,1600. 0*LAL3) 

******************* 

****BUDDEB  PATH**** 

****•*•**•*•**•*•*• 

ESB1=LAL2+BAL2+CSS*2.3 
ESB2=ZHOID(IME40,ESE  1) 

T SA*  ( 1  .0+YK7*  (1 .0-YK8i) 

T5E*  (YK7  +  1.0)  *  (1. 0-YK8) 

T5C=1 . 0- IK8 

8S 53*Z LAG  (BSB2,  Y5A,T£S,  T5C.  IS  E40,BS8  2Z1,3SB3Z1) 
8SB4*asa3«F30  (CC,ES)  *F3  8  <  ALPH  AT,  El,  PS) 

IV2*  F4  2 ( ES84 1 

I V 3D *Z  HO  ID  (I  HP 40  , EG  3)  *ZH CID  (I  E £40,  ?G  3)  * E 1 08  ( QC) 

IV3C  =  TG8»f96  (CC.fS) +HT7- LIMIT  (-30l0, 30.  0,YV33) 

I73B*LIM IT (-30. 5, 20.0.Y73C) 

T73A*ZHOTCH (T73E..13876 . . 27752 . . 1 387 6 .-.79 14 , . 34642, IMP40, 
iV3EZi,Y72BZ2,  Y7  3AZ1.YV3AZ2) 

I73-ZLAG  (T73A,. 44  637, -. 4  4 837,.  10326, IBE40.T7  3FZ1 ,T73Z1) 

8C1=Y72+YV3-Pt6 

RE2*aD  1*F45  (OC.PS) 

HC3*LIMIT  (-30. 0,30.  0.BD2) 

88D1*LIMIT(-3C.O,30.C,SD3) 
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RBD2=A  V4C80  (3 EDI  ,  I  HE  40) 
BBC3=6nT2B  (RATE8U,3HE2) 


* 

* 

* 

* 

« 


RREBEF  =CHFXPL  1 C . 0  ,C .  C,0 . 7 ,4 0 . 0, 1600.  0*EBD3) 
LBC1=LIMIT(-3C.O, *0. C.RD3) 

LRD2=A74C80(LRD1  .IBEuO) 
LBC3=0NTZE(RATE8Q,LHD2i 

Li£DE?  =  Ct!FXPL <0.0,O.C,0.7,4Q.C,1600.  Q*LRD  3) 

»***»*******»**»*#•***«**»*»*♦*« 
*****LEA£ING  EDGE  FLAP  PATH***** 

«****•»»»*»***»»*»*»«*»*»*»»»»»* 

L£2*ZINT  (AL?HAT.C.0625,0. 9375,IBF40,  IE2Z1) 
LE2T=ZH0L£ (IHE20.IE2) 

LE2T1=F27(LE2T,SI) 

LE2T2=>LI2IT(0.  C.?5b  ICC)  ,LE2T1) 

IE3=LI HIT  <0.0,E*9U(Bi),LE2T2) 

L24»?ATLS  (LS3,IHP20. 1B.0*RATE2C) 
LE5*AV2080(LE«.ISP26) 

DIE  1*3  V2*F93  <RI,PS,  NZAF) 

CLE2? 1 =L  25+3 . Q 

DIE2=DEA£S?<-£IE2P  1,  DLE2P  1,  DLE1) 

DLE3*DL2 1+DLS2 
B  t£1=  LE5  *ELE3 

R LE2*LIHIT(- 3. 0,3  3.0. RLE  1) 

HL£3=RATLR  RLE2,IR?80, 18. D*BATE8Q) 
EtEUsLiail (-3 .0,33.0  ,BLE3) 

HIE5=QNTZE  (RAT  £3  0,31*4) 

RLEDEF=CHEX?L  ( 0 . 0 ,0 . 0 ,0 . 7 , 40 .  C.B.LZ5*  16C0.0) 
LtE1=LZ5-DLZ3 

LI£2=LIHIT  (-3. 0,3  2.  O.LL  31) 
LLE3=aATIH(LL22,IHPSC,18.0*RAT£8Q) 

LIE4=  LIHIT (-3. 0,32. O.LL  33) 

"  _  ?  c  /  o  %  rrr/iA 


«0 .  0 , LIES*  1600.0) 


A.  Li  X  .  .  X  I  J  •  W  t  -  -  •  ^  f  -L  J  I 

LLE5=QNTZB  (S  ATEfiO  .Ilia) 

LLEDEF=C2PX?L  <C. 0,0. 0,0.7, 

«*»»****,**«»<»»**»*»»»**»»»»»»** 

•••••TRAILING  EDGE  FLA?  PATH***** 
*«******»•»»**•»*****•***•*»*»*** 

T E2=ZI NT  ( ALPHA!, 0.031 25, 0. 96875, IHP4Q,TE2Z1) 
TS2T=ZH0  ID  (IHP20,TE2) 

TE2T1*F2tt<RI,IE2I) 

T£3*ll SIT  (0.  5,125  (QC)  ,TE2T1) 

T£4=AV208C (TE3  -ISP20) 

0TE1=RATtHfHV2*E711, iHPSO, 18. C*R4TE8C) 
R1E1=TS4-DTS1 

BTE2=LIHIT  (-8. 0,45.0, BT El) 

HIE3=QNTZH  (SATE30  ,BTZ21 

RTEDEF=CRPX?L  (0.0 ,0.0,  ’ . 7 , 40 . C, BTE3*  16CQ.0) 
LIE 1*TS4  *ETE 1 

LTE2=LIHIT(-3. 0,45.0,LTE1) 

LIS3=QNTZR (RAIE80  ,LIE2) 

ITEDEF=C2EXPL  (0.0, 0.0,0. 7, 

TEB  H IN  AL 
HETRCO  HKSFX 

Tia.EB  FISTia*tt  .00,??.DEL*0. 25 ,  OOTDE L* 0. 25,  D  EI‘ 

TITLE  BUECER  3ESBCNSE  TO  */-20LE  TAR 
PHINT  3RECEF , L3DDEF 
PAGE  XYELCT 

OUTPUT  II BE. 3 LEES?  (-5 . 0. 5 . 01  ,  LIEDEF  (- 5 . 0  ,5  .  Ot 
LAEEL  OIFE.L.  2.  FLAP  RESPONSE  10  V-10LB  ROLL, 
LAESL  3*1. 40.H*500CFT 

OUTPUT  TIRE.a TSEEF  (- ?  . 0, 5 . 0)  .  LIEDEF  (- 5 . 0  ,5 . 0) 
LAEEL  DIFF.T. E.  FLAF  RESPONSE  TC  */-10LB  BOLL, 
LABEL  a*1 ,UO,H*500C  FT 

OUTPUT  TIRE.a BCEEF  (-5.0,  5.  0)  .LBEDEFf-  5.0,5 
LAEEL  B  UDDER  REEPCSSS  TO  */-20IE  YAW,  WITH 
LAEEL  a«0.81, H*5000ET 
END 
STOP 

FUNCTION  F4(QC,PS) 

REAL  LIBIT 
F4I1»d.3€-0.0C12*CC 
F4T2«LiaiT(20  0.,  1C0  0.,?S) 

F 413*0. 2-C. 00 03*F4T2 


,40  .  C.LTE3*  1600. 0) 


*0.0125 


WITH 


WITH 


0) 

YG, NY, RG 


1  ~ 


1 


BG  FEEDBACK 
RG  FEEDBACK 
FEEDBACK 
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10 

20 


10 

IS 


10 

20 


F4T4=0.  1-0.00  0  1*F4T2 
I?(.SOT.F4T2.LE.5C0.)  GO  10  10 
1 4T5  *  F4T3 
GO  TO  2C 
CONTINUE 
F4T5*P4T4 
CONTINUE 

F4=LIMIT  (FUTS  ,C.3  ,F4T1) 

RETURN 

END 

FUNCTION  16  (BI.ES.a7  11) 

REAL  LIMIT 

F6T1=0 .22«LIMIT(C.O.  1.  12.  RI1  + (1.0UE-  4)  *FS 
F6T2«LIS IT (0.0,0. «3 1 .F6T1) -0.2S*S711 
F6-.  5  1  1-IIMIT  (£.3  11,0.4  31  ,P6 12) 

RETURN 

END 

FUNCTION  F7(QC,PS) 

REAL  LIMIT 

F7T1=0 .0  5*LIMIT(2S.0, 125.  O.QC)  +3. 75 
F7T2=*QC- 325.0 
F7?3=°S-628. 0 

F7I4=*(LI MIT*(-  325.  ,325.,F712)  * IISIT  (- 336  . ,  0 . ,  P7T3)  )  *  (2.  7 IE- 5)  ♦  1 0 . 0 
IF  (.NOT.  CC.  LE. 225^0)  g6  TC  10 
.7*  F7  T  1 
GO  TO  15 
CONTINUE 
F7*F7T4 
CONTINUE 
RETURN 
ESC 

FUNCTION  FIO(QC) 

SEAL  LIMIT 

F1CT1*1 .C-O.OOC42S*CC 
F 10=* LI  MIT  <0.  3  5.1.C,?10T1) 

RETURN 

ENC 

FUNCTION  F12<ai,PJ) 

REAL  LIMIT 

F12T1*RI**2*9.625-.0  25*  RI ♦ 1 .0 
F 12T2*  PS  *7.  96  9  E-4  *0 .  84 


00396-1.  13) 


F 12T4*  LIMIT  (1.C.8.0.F12T6) 
IFf.NOT.SI.GT .0.5)  GC  TO  5 
F 12  =  E12T4 
GO  TO  6 
CONTINUE 
F  1 2  =  F12T3 
CONTINUE 
RETURN 
ENC 

FUNCTION  FI 3  (RI, PS) 

BEAL  LIMIT 

PS13L*LIKIT(2CC.C,21 16.  0  ,  ?S) 
F13T1=  (4 45.0 ♦LIMIT  (242.  0,  o23.  O.PS 
IE  (.NOT. PS13L.LE. 800.0)  GC  TO  10 
F13T2A*0.054+PE13L*(3. 59E-4) 
F13T2E*-0.  174S-ES13L*(8.4E-5) 
GO  TO  20 
CONTINUE 

F12T2A*0.C0122*PS13l-0.637 
F  13T23=0  .  192-ES13I*  (S.42E-4) 
CONTINUE 

- -  - .34E-5) 


13L) )  *  (2.39E-4) 


3T2  8*81 13L*F13  T2A) 
_  ...  .  _  1  3T  3) 

RETURN 
£  M  C 

FUNCTION  F17  (AIPHAT) 

SEAL  LIMIT 
F17T1*ALPHAT-13.0 
F 17*LI MIT  (0.  0 , 4.  0  ,F  17T1) 


70 


RETURN 

END 

FUNCTION  F22  ( CC) 

f (-SOO- C+LiaiT  <800.  C, 900.0, QC)) 

RETURN 

f  Mf 

F  UNCTION  F23  (ALRH1T) 

F§3*3- iSii-0. 1429*1.1  HIT  (15.0,21.998,  AIEBAT) 

BETURN 

END 

FUNCTION  F24  (SI,  JIPBST) 

?Iul:ii^?38-20.51*LiaiT(0.27,0.66,HI) 

F2  41,2*32. 76-36 ,C*£I3IT70.66,0.91, SI) 
F2uSl=LI!IT(O.C,1COOO.O,  ALPHA!).  „  „„  „  , 

F24T2*  ALFHAf-  ( li .  E  76  9-7 . 6  92  3  *  tl  SIT  ( 0 . 27 , 0 . 9 1  ,BI) ) 
F24I3=-2.C*LIMII  (C.C. 10000.0, **4T2) 

F24T4=  1.  4*  (F24T  1*F^473) 

IF  (.NOT. ai. 31. 0.66)  GO  TO  10 
F24  L=  E24L2 
GO  TO  15 
CONTINUE 

F24L=F24L1 

CONTINUE 

F24*LIMIT  (0.0 ,F24I,F24T4) 

RETURN 

END 

FUNCTION  F25(CC) 

F25=47^636-0.05  106*LIMIT (600.  0,835.0,00 

RETURN 

END 

FUNCTION  F27  (AIrHAT ,RI) 

F27=U  328«(ALPflAT*7.a584-17.8  6»tiaiT  (0.44,0.  63,81)  ) 

RETURN 

END 

FUNCTION  F28(CC) 

j28=44^55  1-0 . 0  40 5 8 ‘LIMIT  (  26  0.  C,  950.  0  ,QC) 

RETURN 

END 

FUNCTION  F29U  (RI) 

F29U=87?Ie25-76.2S*LIMIT (0.7,  1.146,81) 

RETURN 


END 

FUNCTION  S30(CC,?S) 
REAL  LIMIT  , 


F30T1- 1.  0803-  (LIMIT  j  140 . 0,80  0 . 0,gC)  1  *0-0012878 

F  3CT2*  1 . 7  C35-  (II  MIT  (242 . 0  ,4  9  8 . 0 ,  *S)  )  *0.  00  33203 

F30=LIMIT  (F30T2 ,0.9  ,  F30T1  f 

RETURN 

END 

FUNCTION  E31  (Sl.ES) 

l!lTliaiMIT(0.0,2116.0,PS)  ) ‘(6.536E-5)  c  _ 
F  3  1T2=  1. 449- Uiai  ill  455.0  ,21  16.0  ,  PS)  )  *  (5.  IE-4) 
lI(.NOT.  RI.lI.  0.717)  GO  TO  10 
F31T3«F31T2 
GO  TO  20 
CONTINUE 
F31T3*0.  112 
CONTINUE 
7  *  *  C  111 

F31T5*  itlHIT  (-0.717,0.66  2,F3114))*F31T3 
F31T6*g.224  7-?31T1-r 3  ITS 
filiLIMIT  (0.0,0.  1 6.F21T6) 

RETURN 
£  NC 

HSS88.5^SJt"’ 

RETURN 

END 


10 

20 


10 

20 


30 

40 


50 

60 


IT  (0.0,0. 8,  SH  *PS*  (1.69E-4)  -.788 

_  ,  0.0- LIMIT (440.  0, 640. O.rSH  *7.462-3 

F36T3-  {LIMIT  (0.7, 1. 9,31)1  *54  0 .0*955  .  O-QC 


FUNCTION  F34  (  ALPHAT) 

REAL  LIMIT 

F34T1=  1.5-0.  2*  AES  (ALPHAT-2.5) 

F 34-LIMIT  (0.0 ,1.0,?3  4T1) 

PETOaS 
END 

FUNCTION  F35  (  ALFHIT) 

HEAL  LIMIT 

F35T 1-0. 955-0. 03 25* ALPHA! 

F35T2-C. £+0. 04221 *ALPHAT 
IF  (.NOT.  ALPHAT. LE. 6.  C)  GOTO  1C 
F35T3-F35T2 
GO  10  20 
CONTINUE 
F35T3-F35T1 
CONTINUE 

F35-2. C*LIMII  (C.175,0.S,F 35T3 ) 

EETtJBN 

FUNCTION  F36  (RI.CC,  FS1 
REAL  LIMIT 
P36T1-  LIMIT 
P36T2* (44 

F36T3- (LIMIT (0.7, 1.9,31) 1 
F36T4-  ItlMIT] 0.0, 200.0,  F36T3)  )  *0.005 
F36T5=F3€T1*,F36l£*1.C 
F 3616=  LIMIT  (0.0, 1.0,  F36T5) 

F36-AMIN 1  (F3  6T4.7  36T6) 

EETURN 

END 

FULCTION  {37  (HZ 2) 

REAL  LIMIT 

F37=2.  S-0.5*LIBIT  (3. 0,5.  Q  ,HZA) 

RETURN 

£  N  D  • 

FUNCTION  F38(ALPHAT,RI,?S) 

HEAL  LIMIT 

FSRI=LIM IT (0.225,0.549,211 
IF  (.NOT.  ALPHAT.LT.  16. 0)  GOTO  10 
FS— 0.243 C  8 *F SSI* 0.03 967 
AOAM—10.  0 
A0A3  —  0.45 
GO  TO  20 
CONTINUE 

FS—0.  13846*FSEI-0.0  1774 
ACAM—  5. 467 
ACAE=-0.  025 
CONTINOE 

F38T1- LIMIT  (-5. 0, 20.4,  ALPHA!) 

F36T2-LIMIT(-5.C,25 .O.ALPHATL 
IF(.  NOT.  SI.  LE.  0.  354)  GO  TO  30 
F38T3-F38I2 
LL= -0 . 98 
GO  TO  40 
CONTINUE 
P36T3-P36T1 
II— 1.  2 
CONTINUE 

F38T4-  (F38T3  +  ACAM) *FS*AOAE 
F  38TS-LI MIT(LL.0.3.F3  8T4) 

F38T6-  0.048-0.  04 9 6 -LIMIT  (-5. 0 ,25 .0,  A LPH AT) 

F3  8T7- LI  MIT  (-5, 0,25.0,  ALP  HAT)  *0 . 0  1*LI3IT  <240 . 3 ,390.  ( 
F38T9-LI MItIo .C,2.5,F38T7)-2.5 
F38T9*  BI-0.  72 

F38T10-0.C168*LI3IT(C.Q,2.08,F38T9) 
F38T11-F28T10*F38T8*F38T6 
IF  (.NOT.  A  LPH  AT  .IS  .  1 0 .0)  GC  TO  50 
F38-F3ET1 1 
GO  TO  60 
CONTINUE 
F38-F38T5 
CONTINUE 
SFTUSN 
END 

FUSCTION  F39 ( AIPHAT) 

SEAL  LIMIT 


,?S)  -3.9 


72 


F 39* (LIMIT ( 1 3. C, 25. C, ALPHA?)  )  *0.0833  3-1.08329 

BETUBN 

END 

FUNCTION  E40  (RI.BS.QC.PIQ) 

P4 011=  3? i?-3. 0*11  MIT  [0.05, 0.83, 3 XI 
740T2=C.65625-.0C13  125*1.1  MIT  (500. 0.  1  800 . 0  ,  PS) 

F40T3* -0.26 177*  (<3.£3  5E-u)  ‘LIMIT  (5G0 .  C,  190  5.0,??) 

F40T«=  II SIT (0.0.  1C5.  C,QC-  335.C1 

F40T5*  (F4CT4*  (-EIC*  (1.17  4  28  E-61*  (1.5  238  E-3)  )  )  *0 . 475 -  (6 . 5E-  4)  *PIQ 
F4  0T6=  F40T3*  (E4QT2*LIMIT  (  0.  7  5.0 .85,  3  I)  ) 

F4CT7=F4CI4*]S  .67  46  E -4*  1 .  98  4  1  f-7*EIQ)  ♦  (0 .  1 69  23-3 . 84 6 1 5  E-5 *  P IQ) 
F4CT8=  <0.1692  3-3. 864  152-5  *P  IQ)  ♦  (F40T  4*  ( 1 . 672472-3 -9 . 29  152  E- 7*FI3)  | 
IF  (.SCT.BI.GE.0.75)  GO  TO  20 
F40T9=F40T6*F4CI1 
GO  TO  30 

20  CONTIS  OE 

P40T9*E40T8 
30  CONTIS  US 

IF  (.NOT.  (QC.  GE.  440.0.  AND.  SI.  GE.  0.751)  GO  TO  40 
F40T1 C=F40I6 
GO  TO  50 

40  CONTINUE 

F40T  10=F40T7 
50  CONTINUE 

If  (.  NOT.  Pig.  GT. 980.0)  GO  TO  6C 
FU0  =  F  UCT9 
GO  TO  80 

60  CONTIN  UE 

IF  (.MCT.Eig.LE.500.C1  GO  TO  7C 
F40=F  4CT5 
GO  TO  80 

70  CONTINUE  „ 

F40  =  F40T  1 0 
80  CONTINUE 

BETUBN 
END 

FUNCTION  F41  (ALPS AT, 3X1 
REAL  LIMIT 

F41T1*  (LIMIT  (8. 0,22. 0.ALPHAT1  -6.01  *0  .05  3571 

F«  1T2=i.  2353  3*LIMIT(.  27..  48.311  -.009  971  *  ( LIMIT  (8  . ,  2  2.  ,  ALPH  AT ) -3 . ) 
F4  1T4=0.  5875- (LIMIT  <0.2 5.0.59, Si)  )*1 .25 
IF  (.NOT.  AIPHAT.L2.8.  0)  GO  TO  10 
F41T3»1.0-F41T1 
GO  TO  20 
10  CONTINUE 

F4  1T3=  1.  0-F41T2 
20  CONTINUE 

F4  1=LI SIT  (F4  1T4 , 1.0,  F41T3) 

HETUHN 

ENC 

FUNCTION  F«2(IN) 

BEAL  LIMIT,  IN 

P 42*  (A as  (LIMIT  <- 18. 3 3. 18.  33  ,IS))  *0.0479*0.76  1) 

*  ‘LIMIT  (-18.  i3,  16.  33,  IN) 

BETU3N 

ENC 

FUNCTION  F45(CC,P£) 

SEAL  LIMIT 

F45T1=(LIMIT(5C0.C,10g0.0,gC) -500.01  ‘O^OI 

PU5T2*  (LIMIT  (625. Oj  1450. 0,?S  )  *0.  00  1  2 1  ‘  1-0 . 75758 

F  45*1 . C-?45T1*F45I2 

BETUBN 

ENC 

FUNCTION  F68(QC) 

BEAL  LIMIT 

F66*-O.OC2977*  (-460 . 0* LI  MIT  (2  60. 0,4  8  0. 0  , CC) ) 

BETUBN 

END 

FUNCTION  F90  (BI,  PS) 

H^i=ili?o-ps 

P90T2-RI-C.72 

F90  =  (LIMIT  (0. 0,386.  0.F9  0T11  ‘LIMIT  (0.  0,2.08,F90T2)  )  *0.0  2055  1*16.5 

BETUBN 

ENC 


73 


10 

20 


10 

20 


FUNCTION  F93  (FI.FS.NZAF) 

BEAL  L IMIT,N  ZA  F 

F93T1* (RI*0.00018»PS-0.69) *0. 4  2857 
F92T2=PS-628 . 0 

F93T3=  1-0-0.  1*LIMIT(G. 0,5.0,  NZAF) 

F93T4=LIMITf 0.0,0. 06 ,F93T1) 

F9  3X5  =  (6.3694E-3)  »LI».IT  (0.0, 157.0 ,F9  3T2) 

F93=F9  3T4*F9  3TS*  r  S  3T  2 

RETURN 

END 

FUNCTION  F96  <QC,P£) 

BEAI  LIMIT 

F96T1*1.  S«-0.0 C44UIJ ‘IT  (0.0, 1  68.0, QC) 

F96T2=CC*0. 00088*0. 36 
IF  (.NOT. CC.LE  .500.0)  GO  TC  10 
F96  13«F96T1 
GO  TO  20 
CONTINUE 
F96T3=F96T2 
CONTIN  UE 

F96T«  =  0.  000777  2*  (6  28.0-LIMIT  (  242.  0, 6  28.  0 ,  ?  S)  ) 

F96T5=P96I3*F96I4 

F96*LIMIT  (0.0,  1.  1.F96T5) 

RETURN 

ESC 

FUNCTION  F101  (BI.FS.NZ1F) 

BEAL  LIMIT. NZAF 

F101T1 =R I-Q . 5  6 

F101T2-LIKIT  (O.C,5.0,F101T1) 

F1C1T3*P£-498.0 

F101T4=F  10  IT  2*  LIMIT  (0.0, 835.  O.F101T3) 

F101TS=79E.54-C.35S48*PS 

F 1 0 1T6 “LIMIT  (  165. 2,275.  5,  F101T5) 

F1C1T7=(5.4u46E-5|  *LIMIT(0.  0  ,  F  10  1T6,  F  10  1T4) 
F1C1T8=AE£<NZAF) 

F 10 1T9» LIMIT  (0.0,5.0,F10  1T8) 

F1C1=?101T7*F1C1T9 
RETURN  ' 

ENO 

FUNCTION  F107(CCI 
REAL  LIMIT 

F  10711=  /  (-5.  714E-7)  *ABS  (QC-7S0.0)  )  ♦  (8.4E-4) 
F1C7=LIMIT(O.C,10COO.O,F107T1) 

EETURN 

END 

FUNCTION  F108(CC) 

REAL  LIMIT 

F109T1=LIKIT  (5  0.0,390.0,QC) 

F 1 08=  (F  1  0811**2.0)  *7.48 E- 8+0.  01557-F  108T1*6.  4352-5 

RETURN 

END 

FUNCTION  F 1 12  (El) 

REAL  LIMIT 

IF  (81.  IT.  0.  28.CR.FI.  GT.  0.76)  GC  TO  10 
F 1 1 2  =  4 .0 
GO  TO  20 
CONTINUE 
F1 12=0 . 0 
CONTINUE 
EETURN 
END 

FUNCTION  F1 1  3  ( 1LPEAT) 

REAL  LIMIT 

F113=0  .  16  66  7*  LIMIT  (12.0,  18.0,AL?HA7)  -2.0 

RETURN 

INC 

FUNCTION  F 1 1 4  (RI) 

REAL  LIMIT 

F1  14=0. 89286*  LIMIT  (0.  16 , 0.30  ,  FI)  -0.  1  4  28  6 

EETURN 

END 


END  JOB 
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